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Abstract 

Two-jet production in deep-inelastic electron-proton scattering has been studied for 
160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 and 0.04 < y < 0.95 with the ZEUS detector at 
HERA. The kinematic properties of the jets and the jet production rates are presented. 
The partonic scaling variables of the two-jet system and the rate of two-jet production 
are compared to perturbative next-to-leading order QCD calculations. 
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1 Introduction 



Deep-inelastic neutral current scattering is described by the exchange of a virtual boson (7*, Z°) 
between the electron and a parton in the proton. In the naive Quark-Parton-Model (QPM), 
this process leads to a 1+1 parton configuration in the final state which consists of the struck 
quark and the proton remnant, denoted by "+1". Higher-order QCD processes contribute 
significantly to the ep cross section at HERA energies: to 0(a B ) these are QCD-Compton 
scattering (QCDC), where a gluon is radiated by the scattered quark and Boson-Gluon- Fusion 
(BGF), where the virtual boson and a gluon fuse to form a quark-antiquark pair. Both processes 
have 2+1 partons in the final state, as shown in Fig. [T]. 

Jet production in Deep-Inelastic Scattering (DIS) has been studied by a fixed target experiment 
(E665) at a centre of mass energy, y/s, of 31 GeV and at negative squared momentum transfers, 
Q 2 , of order 10 GeV 2 Jl||. At the much larger centre of mass energy of 296 GeV at HERA, 
jet structures are much more visible 0, E|. It is therefore possible in this energy regime to 
determine the Q 2 dependence of the strong coupling constant, a 8 , in a single experiment by 
measuring the rate of two-jetQ production ||. 

Such measurements require a detailed understanding of jet properties. In this paper we study 
whether the two-jet system has large enough invariant mass and whether the jets have suffi- 
ciently large transverse momenta for perturbative calculations to be applicable. Using Monte 
Carlo simulations of parton showering and hadronisation to correct for higher-order and non- 
perturbative effects, the underlying parton dynamics and the jet rates are compared to next- 
to-leading order (NLO) calculations. 

We study jet production in the kinematic range 160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 and 
0.04 < y < 0.95, at an average squared hadronic invariant mass (W 2 ) of 14000 GeV 2 , where 
x is the "Bjorken x" variable and y denotes the energy fraction transferred from the incoming 
electron to the proton in its rest system. The choice of this kinematic range was based on the 
following expectations: at high Q 2 , jet structures should be more pronounced and hadronisation 
uncertainties should be reduced. Furthermore, in the region x > 0.01, the phase space for jet 
production increases. Also, at x > 0.01, the theoretical uncertainties in the jet rates due to 
different parameterisations for the parton densities of the proton are small ||. 

The data were collected with the ZEUS detector in 1993 and correspond to an integrated 
luminosity of 0.55 pb _1 . A brief description of the ZEUS detector is given in section |2|. The 
event selection is explained in section [| In section ^ we describe the Monte Carlo simulation 
and the theoretical calculations. The JADE algorithm 0] which is used to relate parton- and 
hadron-level jets in this paper is described in section [| The kinematic variables for two-jet 
production are defined in section [| In section [7| we study the hadronic energy flows, the 
pseudorapidity distribution of the jets in the detector, the partonic scaling variables and the 
jet production rates. Conclusions are given in section |8|. 



2 The ZEUS detector 

The experiment was performed at the electron-proton collider HERA using the ZEUS detector. 
During 1993 HERA operated with bunches of electrons of energy E e = 26.7 GeV colliding with 
bunches of protons of energy E p = 820 GeV, with a time between bunch crossings of 96 ns. 

1 In this paper we refer to 2+1 jet production as two-jet production. 
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ZEUS is a hermetic multipurpose magnetic detector that has been described elsewhere ]7|, 
Only components relevant to this analysis are mentioned here. 

The hadronic final state and the scattered electron are measured by the uranium-scintillator 
calorimeter (CAL). It consists of three parts, the Forward (FCAL), the Rear (RCAL) and the 
Barrel Calorimeter (BCAL)^. Each part is subdivided longitudinally into one electromagnetic 
section (EMC) and one hadronic section (HAC) for the RCAL or two HAC sections for BCAL 
and FCAL. Holes of 20 x 20 cm 2 at the centre of FCAL and RCAL accommodate the HERA 
beam pipe. In the XY plane around the FCAL beam pipe, the HAC section is segmented in 
20 x 20 cm 2 cells and the EMC section in 5 x 20 cm 2 cells. In total, the calorimeter consists 
of approximately 6000 cells. In terms of pseudorapidity r\ = — In tan |, the FCAL covers 
4.3 > T) > 1.1, the BCAL 1.1 > 77 > -0.75 and the RCAL -0.75 > 77 > -3.8, assuming the 
nominal interaction point (IP) at X = Y = Z = 0. The CAL energy resolution as measured 
under test beam conditions is ue/E = 0.18 /VE (E in GeV) for electrons and ue/E = 0.35/ y/E 
for hadrons. The timing resolution of the calorimeter is less than 1 ns for energy deposits greater 
than 4.5 GeV. 

The beam monitor scintillation counter (C5) was used to measure the timing of the proton and 
electron bunches. The event vertices were determined by drift chambers surrounding the beam 
pipe: the Vertex Detector (VXD) and the Central Tracking Detector (CTD). The resolution of 
the Z coordinate of the primary vertex is 4 mm. 



3 Event selection 
3.1 General selection 

During the 1993 data-taking period about 10 6 events from the DIS trigger branch were recorded. 
A more detailed description of the DIS trigger conditions and some aspects of the event selection 
can be found in |J. The trigger acceptance was essentially independent of the DIS hadronic 
final state with an acceptance greater than 97 % for Q 2 > 10 GeV 2 . In order to select DIS 
events, the following set of cuts was applied: 

• A timing cut required that the event times measured by the FCAL and the RCAL were 
consistent with an interaction inside the detector. This cut strongly reduced beam-gas 
background. 

• A scattered electron candidate had to be reconstructed with an energy E e i greater than 
10 GeV to ensure high purity of the electron sample. 

• The position (X, Y) of the scattered electron in the RCAL had to lie outside a square of 
32 x 32 cm 2 centred on the beam axis, to ensure the electron was fully contained within 
the detector and its position could be reconstructed with sufficient accuracy. 

2 The ZEUS coordinate system is defined as right handed with the Z axis pointing in the proton beam 
direction, hereafter referred to as "forward" . The X axis points towards the centre of HERA, the Y axis points 
upward. The polar angle 9 is taken with respect to the Z direction. 
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• Events were selected by the requirement 35 < 5 = J2i Ei(l — cos^) < 60 GeV, where 
Ei, 9i are the energy and polar angle (with respect to the nominal IP) of the calorimeter 
cells i. For fully contained events 5 ~ 2E e = 53.4 GeV, where E e is the electron beam 
energy This cut was applied in order to remove background due to photoproduction and 
beam-gas interactions. 

• The Z position of the event vertex was reconstructed from the tracking data. Events were 
accepted if the Z position was inside ±75 cm of the nominal IP. 

• Events from beam halo muons, cosmic rays and QED Compton processes were identified 
and rejected by suitable algorithms. 



3.2 Kinematics 

Because the ZEUS detector is almost hermetic, the kinematic variables x, y and Q 2 can be 
reconstructed in a variety of ways using combinations of electron and hadronic system energies 



and angles [10 



1. The electron method in which the kinematic variables are reconstructed from the energy 
E e i and angle 6 e i of the scattered electron. 

2. The Double Angle (DA) method in which the angles of the scattered electron (9 e >) and 
of the hadronic system (j H ) are used. This method reduces the sensitivity to energy 
scale uncertainties. The angle jh corresponds to that of a massless object balancing the 
momentum vector of the electron to satisfy four-momentum conservation. In the naive 
QPM 7# is the scattering angle of the struck quark. It is determined from the hadronic 
energy flow measured in the detector using the equation 

_ (J2 h Px) 2 + (E h py) 2 -(J: h (E- P z)) 2 m 
C ° SlH (E h Px) 2 + (E h PY) 2 + (E h (E-p z )) 2 1 1 

where the sums, Ylhi run over a ^ calorimeter cells h which are not assigned to the scattered 
electron, and (px,Py,Pz) is the momentum vector assigned to each cell of energy E. The 
cell angles are calculated from the geometric centre of the cell and the vertex position of 
the event. 



The Jacquet-Blondel (JB) method [11| in which the kinematic variables are calculated 
from the reconstructed hadronic final state using the momentum vector (px,PY,Pz) con- 
structed from the energies E and angles 9 of the calorimeter cells which are not assigned 
to the scattered electron. 



The electron method gives better resolution in x at low Q 2 while the DA method is less sensitive 
to the calorimeter energy scale and gives the better mean resolution over the whole x-Q 2 plane. 
The JB method is used to calculate the visible hadronic energy in the events which will be used 
in the jet definition. Complete formulae for calculating the variables x, y and Q 2 are to be 



found in fp7Q]| . When it is necessary to distinguish which method has been used the subscripts 



'DA 1 or l JB J will be used on the variable concerned. 



3 



The data sample used for this analysis had to satisfy the following cuts: 

• In the DA method, in order that the hadronic system be well measured, it is necessary 
to require a minimum of hadronic activity in the CAL away from the beampipe. For this 
purpose the value of yj^ = J2hE(l — cos9)/2E e had to be greater than 0.04. 

• The value of y e = 1 — — cos9 e i) had to be less than 0.95. This cut rejects photo- 
production background. 

The kinematic variables Q 2 and x were determined with the DA method. Starting from 
28200 events with Q^ A > 10 GeV 2 , 1020 events remain after selecting the kinematic region 
160 < Qda < 1280 GeV 2 and 0.01 < a; da < 0.1 for this analysis. The remaining photoproduc- 
tion background in this high (x, Q 2 ) data sample is <,1%. 

The corrected distributions are corrected for detector and acceptance effects and given in the 
kinematic ranges 160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 and 0.04 < y < 0.95. To investigate 
the evolution of jet structures with Q 2 , a low (x, Q 2 ) data sample of 4230 events is used with 
10 < Q 2 < 20 GeV 2 , 0.0012 < x < 0.0024 and 0.04 < y < 0.95. 



4 QCD calculation and event simulation 
4.1 LO simulations 

The detector effects were simulated using events generated with the LEPTO 6.1 Monte Carlo 



program |12| based on the electroweak scattering cross section. We have studied the process 



ep — ► e'X considering 7* and Z° exchange. 

The leading order (LO) matrix elements (ME) are used to simulate the QPM process, 7*+g — > q, 
the QCDC process, 7* + q — > q + g, and the BGF process, 7* + g — > q + q (Fig. [[J). In order 
to obtain finite cross sections for parton emissions, a minimum invariant mass, m 2 j = y m i n W 2 , 
is required for all pairs of final state partons i and j. Higher order parton emissions, which 
are calculated in the leading log approximation (LLA) of perturbative QCD, are simulated 
by the Parton Shower (PS) model. The struck parton can radiate partons either before or 
after the interaction. The amount and the hardness of the radiation depends on the virtuality 
(mass) of the partons. In order to simulate the hard emission of partons and the higher-order 
parton showers, a combined option (MEPS) exists. The value y m i n = 0.015 was used. The 
parameterisation of the parton distribution functions was the MRSD'_ set which has been 
shown to describe reasonably the HERA measurements of the proton structure function, F2 0. 

The event generation included the effects of initial- and final-state photon radiation which 



were calculated with the program HERACLES 4.4 [[Uj]. The simulation of the detector used a 
program based on GEANT 3.13 H, ITJ. 
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4.2 NLO calculations 



The LEPTO 6.1 Monte Carlo event generator uses the exact 0(a s ) matrix element (ME) 
and the parton shower (PS) in the leading log approximation. It does not include the NLO 
matrix element calculation. However, the NLO corrections to the 2+1 jet cross section due to 
unresolved 3+1 jet events and due to virtual corrections are significant |J. These corrections 
are included in the program DISJET of Brodkorb and Mirkes [I5| and in a similiar program by 



Graudenz called PRO JET 117 



Both programs calculate cross sections at the partonic level in NLO as a function of x and Q 2 . 
PROJET 3.6, in addition, provides cross sections in terms of the parton variables (see section 
^|), but it does not contain the NLO corrections to the longitudinal cross sectionP]. Only the 
exchange of a virtual photon (7*) is considered in the calculations, however the contribution 
from Z° exchange is small (expected to be less than 2% for Q 2 around 1000 GeV 2 11). 



5 The jet finding algorithm 

A jet finding algorithm is necessary to relate the hadronic final states measured in the de- 
tector to hard partonic processes. Results on two-jet production in DIS have been pub- 
lished previously by the ZEUS collaboration based on a cone algorithm using a cone of ra- 



dius R = y(Ar/) 2 + (A0) 2 = 1 §; here the cone variables are the pseudorapidity 77 and the 
azimuthal angle 0. 

In this paper we use the JADE algorithm ||, since it is currently the only algorithm which allows 
comparison to the NLO calculations. The performance of the JADE algorithm in reconstructing 
jets has been compared to other algorithms in ||18|| . The JADE algorithm is a cluster algorithm 
based on the scaled invariant mass 

JADE = 2EjEj(l - COS Ojj) 

for any two objects i and j assuming that these objects are massless. W 2 is the squared invariant 
mass of the hadronic final state and 9y is the angle between the two objects of energies E$ and 
Ej. The minimum y^ of all possible combinations is found. If the value of this minimum 
is less than the cut-off parameter y cut , the two objects i and j are merged into a new object 
by adding their four-momenta and the process is repeated until all yij > y cnt . The surviving 
objects are called jets. 

The JADE algorithm is applied at the parton, hadron and detector levels in the HERA labo- 
ratory frame. At the detector level, the calorimeter cell energies and positions serve as inputs 
to the JADE algorithm. For this analysis, the scale parameter used in the JADE algorithm, 
W 2 , is taken at the detector level from just those calorimeter cells associated with the hadronic 
system, W 2 B — s(l — xj B )y JB . By using Wjb, the detector acceptance corrections for y^ largely 
cancel. 

We have found that the JADE scheme as it is described here leads to smaller hadronisation 



corrections than the Lorentz invariant i?-scheme |19| which uses the exact invariant mass 



for massive objects. For this analysis, the JADE scheme is used in the HERA laboratory frame. 
3 These corrections are contained in a newer version (PROJET 4.1.1). 
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The losses in the forward beam pipe of the ZEUS detector are taken into account by adding 
a fictitious cluster (called pseudo-particle) in the forward direction to which the missing lon- 
gitudinal momentum for each event is assigned. The pseudo-particle is treated like any other 
particle in the JADE clustering scheme. No pseudo-particle procedure is used for the parton 
and hadron levels of the Monte Carlo generator. 



6 Two- jet kinematics 

The differential cross section for two-jet production depends on five independent kinematic 
variables ||20|| , here taken as x,Q 2 ,x p , z and 0*. The variable <fi* is the azimuth between the 
outgoing parton plane and the lepton scattering plane in the 7*-parton centre of mass system 
(CMS) and the variables Lorentz invariant partonic scaling variables. 



The event variable x p is determined by: 

Q 2 _ x 

2p-q £ 



x v 



(0 <x p < 1), 



where £ is the fraction of the proton's (longitudinal) momentum P carried by the incoming 
parton of momentum p (p = £P) and q is the momentum of the exchanged virtual boson 
(Q 2 = — q 2 ). Alternatively, one can write: 

Q 2 

X, 



up 



Q 2 + s 



for massless jets, where y/1 is the invariant mass of the two-jet system. This expression is used 
to determine x p experimentally. The range of values of x p is given by x as the lower limit and 
is a function of y cut at the upper limit: 

Q 2 Q 2 
X Q 2 + W 2 ~ Xp - Q 2 + y cut W 2 - 

The jet scaling variable z is related to the angular distribution of the jets in the 7*-parton 
CMS: 

*=^ = 1(1 -008^) (0<Z<1), 

where 6? et is the polar angle of the jet in the 7*-parton CMS. z is measured from the jet, which 
is assumed to be massless, and reconstructed in the HERA system: 

£jet(l - COS f jet) 



J2iEi(l- COS 9i) 



where the sum runs over the two reconstructed (massless) jets. In the JADE algorithm, the 
total hadronic energy is contained in the 2+1 jet system. Since the remnant jet ('+1') has no 
transverse momentum, Remnant = and therefore Z1+Z2 = 1. The minimum value is determined 
by the value of y cut and is given by: 



_ (1 - x)x p 

Zmin ^ ^ ycut 
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In this analysis, we integrate over fa and study the dependence of the two-jet production on 
x p and z. The transverse momentum Pt of the jets with respect to the 7* direction in the 
7*-parton CMS, can be derived from x p , z and Q 2 : 



lQ 2 -(l-x p )(l-z) 

Xp 



At the 0(a s ) tree level, the singularities in the two-jet cross section are given by plf : 

v 2 1 n ,2\1L2 1 n „2M 11 ^2^2 

7 J/ ' 

'2+1 



1 + £p£ g 



^ (1-X p )(l-^)' 

where z g is labelled specifically for the quark jet in the QCDC process. The x p singularity 
in the QCDC process results in small jet-jet invariant masses being preferred. Both processes 
have a singularity at z — or 1: in the BGF process, it is related to the collinear emission 
of the two quarks and, in the QCDC process, to the collinear or soft emission of the gluon. 
The kinematic properties of the JADE jets and the relationship to these 0(a s ) singularities are 
studied in section ITT 



7 Results 

7.1 Transverse energy flow in jets 

The jet properties were studied at a fixed y cut of 0.02. This value is chosen as a compromise 
between the increase of the higher-order corrections at lower y cut values and the loss of statistics 
at higher y cut values for two-jet events. This choice resulted in 237 events with 2+1 jets from 
a total number of 1020 events in the high (x, Q 2 ) region. The distributions in the following 
sections are normalised to the total number of events, N ev . 

The transverse energy (Et) flows relative to the Z axis in the laboratory frame 

1 dE T , 1 dE T 
; r and — — r 

N ev d(Arn) N ev d(Afa) 

were calculated with respect to the jet directions in the laboratory frame by defining 

Arji = r]i - ?7 jet and A fa = fa - fa ct 

for every calorimeter cell i of 1+1 jet events (Figs. |2]a,b and |^a,b) and separately relative to 
the axes of the higher and the lower t] jets of 2+1 jet events (Figs. |2]c-f and |3|c-f). 

In Fig. the Et flows are compared in two different kinematic regions at low Q 2 (10 < Q 2 < 20 
GeV 2 , 0.0012 < x < 0.0024) and at high Q 2 (160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1). 
The average Et flow increases with Q 2 and the jets are more collimated at higher Q 2 . From 
kinematics (Et) — Q for the electron and therefore the same is expected for the balancing 1+1 
jet. For the two kinematic regions shown here, this corresponds to approximately a factor 5 
increase in (Et)- The influence of the proton remnant is visible as a tail at At] > 1 in the 1+1 
jet events (Fig. ^|b), whereas the tail in Fig. is mainly assigned to the second jet and not 
the proton remnant. 
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Fig. [3] shows the Et flows for the high (x, Q 2 ) region compared to two Monte Carlo data sets 
after the detector simulation, generated with the ME and the MEPS options. The Et flow in 
terms of the energy scale and the shape around the jet direction is well described, especially 
when using the MEPS option. In the ME model, the absence of parton showering results in 
the Et flow being more concentrated around the jet direction. 

Every cell contributing to the Et flow is assigned to one of the jets by the JADE algorithm. 
The Et flow contribution of the jet defining Arj = or A<p = is shown as the shaded 
histogram in Fig. ^. The JADE algorithm is observed to typically assign cells to a jet inside 
the range |Ar/|^,l, however, the algorithm also assigns cells beyond this range for the higher rj 
jet (Fig.|e-f). 

7.2 Properties of two-jet events 

The two-jet properties are studied in this section for the high (x, Q 2 ) region. In Fig. |], the 
distribution of the pseudorapidity r/j et of the two jets is shown. The jets are ordered in rj. The 
higher rj jet is usually found very close to the forward beam pipe. About half of the jet axes 
are reconstructed in cells within 30 cm of the FCAL region near the beam pipe (8<8°). In 
this forward region, the results depend on the description of the initial state parton shower and 
of the target fragmentation in the Monte Carlo generator, as well as on the simulation of the 
response of the calorimeter around the beam pipe. 

Fig. [|a shows that the predictions of the 77j et distribution by the ME and the MEPS models 
describe the data fairly well except for the very forward angles (r/j et > 3.6) where the predictions 
are below the data for both models. In the region 3.6 > r]j Ct > 2.8 the data are better described 
by the MEPS model. Most of the lower rj jets are also found at positive pseudorapidities in 
the BCAL or in the FCAL. The Monte Carlo model predictions for the r]- ]et distribution of the 
lower 7] jet, which is usually well separated from the beam pipe region, gives in general a good 
description of the data (Fig. |]b). 

The distributions of the differences in azimuthal angle, A0j et , (Fig. |5]a) and pseudorapidity, 
A?7j Ct , (Fig. |5|b) between the two jets show that the jets found by the clustering scheme of the 
JADE algorithm are reasonably well separated in rj. Since <fi is calculated in the laboratory 
frame and the electron acquires more transverse momentum at larger values of Q 2 , the two 
jets are not back-to-back in the laboratory frame for the high (x, Q 2 ) data. This differs from 
low-Q 2 jet production, where the jets are typically back-to-back in <p. The data are generally 
well described by the Monte Carlo predictions for A?7j et and A0j et . 

The invariant mass v s of the two jets is shown in Fig. ^|a. The average invariant mass (VI) of 
the jets is about 23 GeV. The yl distribution is reasonably well described by the two models. 
The Pt of the jets in the 7*-parton CMS, determined from the expression in section 6, is shown 
in Fig. |5|b. (Only the Pt distribution of one of the two jets is shown, since \P T 1 |=|Pj < |.) 
The average value {Pt) is about 7 GeV/c. This is sufficiently large to ensure the validity of 
a perturbative QCD calculation. At low Pt, however, the data lie above the Monte Carlo 
prediction. 
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7.3 Partonic scaling variables 



The distribution of z versus x v for the uncorrected data is presented in Fig. |7|. The area 
defined by the curve specifies the kinematic limit for y cut = 0.02 in the JADE algorithm and 
0.01 < x < 0.1. The upper and lower limits on x p depend on x. The upper limit varies in the 
range 0.34 < x p < 0.84 for 0.01 < x < 0.1, respectively. The data typically lie close to the z 
limit, given by z min ~ y cut = 1 — z max , which is close to the singular regions of the cross section 
discussed in section 6. 

The measured x p distribution for the two-jet events is shown in Fig. |8|a at the detector level. 
The bin size has been chosen to reflect the experimental resolution of around 10% in x p . In 
the considered (x, Q 2 ) region, there exists almost no phase space limitation at the upper end 
of the x p distribution, with x p extending up to approximately 0.8 for the high Q 2 data. This 
is however sufficiently far from the x p singularity of the QCDC process discussed in section [(]. 
The average x p is about 0.5, i. e. (s) ~ (Q 2 )- The x p distribution is well described by the MEPS 
and the ME models. 

The uncorrected x p distribution for the low Q 2 data which is shown for comparison shows a 
very different behaviour. It is peaked at small values of x p , i. e. (s) > (Q 2 ), due to the y cu t 
requirement. At these lower Q 2 values, the largest scale is therefore s. 

In the uncorrected z distribution in Fig. |8|b, only the z distribution for the jet with smaller z 
is shown, since z\ + Z2 = 1. At the detector level, the distribution is reasonably well described 
by both the ME and the MEPS model, but a discrepancy exists for z < 0.1 which corresponds 
to jets close to the forward beam pipe. 

The corrected x p and z distributions are compared to the PROJET NLO calculation in Figs, ^c 
and d. These distributions were corrected for detector and hadronisation effects with the MEPS 



simulation. The MRSD'_ parton density parameterisations |y| are used in the calculations 
of PROJET. Apart from the parameterisation of the parton densities, the QCD calculations 
contain only one free parameter, the strong coupling constant a s . For the calculations in Fig. |] 
a value of = 250 MeV was chosen which corresponds to a s (M|) = 0.120, as measured from 



hadronic event shapes, energy correlations and jet rates in e e annihilations [IS 



The shape of the corrected x p distribution and the PROJET calculation is in good agreement, 
however, the rate is underestimated. The corrected z distribution is well-described by the cal- 
culation for z > 0.1, however, an excess is observed in data relative to the PROJET calculation 
in the lowest z bin. This is, however, the region with the largest systematic uncertainties. 
The significance of this deviation cannot currently be estimated due to the experimental and 
theoretical problems associated with the jets in the forward region. The z distributions are 
almost identical for the ME(LO), not shown, and the PROJET(NLO) calculation. The cross 
section rises as z — > 0, because of the low cut-off with respect to the z singularity in the JADE 
algorithm. 

It should be noted that the jets at the MEPS parton level acquire mass in the JADE algorithm 
through multiple combinations of partons from the parton shower (PS) whereas the PROJET 
and the ME jets are defined by massless partons. This effect shifts both the x p and the z 
distribution as defined in MEPS, if x p and z are not calculated in the approximation of massless 
parton jets. 
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It should also be noted that the kinematic cut-offs in terms of the partonic scaling variables 
strongly depend on the choice of jet algorithm. In comparison to the JADE algorithm the K± 
algorithm |22[ leads to a much less peaked z distribution. The K± algorithm was not used for 
this analysis, because the NLO calculations for DIS are based only on the JADE scheme. 



7.4 Jet rates 

The two-jet rate -R2+1 is defined by the ratio 



R 



N. 



2+1 



2+1 



N 2+1 + N : 



1+1 



where N 2+ i and iV 1+1 are the number of 2+1 or 1+1 jet events and R\ + i = 1 — R 2 +i- This 
definition differs from the "usual" definition where the denominator is the total cross section 
including 3+1 and higher-order contributions. It is used to reduce the dependence on the 3+1 
jet rate, which is only calculated at the tree level in PROJET and DISJET. Experimentally, 
the acceptance correction factors for the 3+1 jet rate are large. For y cnt = 0.02 about 4% of 
all events have 3+1 reconstructed jets. Almost no events with 0+1 jets are found in this high 
(x, Q 2 ) interval whereas they are often found in the lower (x, Q 2 ) intervals. 

In Fig. the corrected jet production rates are shown as a function of the jet resolution 
parameter y cut . The correction to the parton level is done using the MEPS model. The 2+1 jet 
rate increases with finer jet resolution (smaller y cut ). The measured jet rates have been corrected 
using a bin-by-bin correction method, where the correction factors are up to 20% depending 
on y CU f At lower y cut values (y cnt — 0.02) the resolution on y^ is approximately 0.01. 

For comparison a full NLO calculation performed with the DISJET program and the NLO cal- 
culation of the PROJET program are also shown. The difference between the two calculations 
is attributed to the missing NLO corrections to the longitudinal cross section in PROJET 3.6. 
The PROJET and DISJET curves are calculated using A® = 250 MeV and the MRSD'_ par- 



ton density parameterisations, as discussed in section 7.3. Choosing other currently available 



parameterisations G3| leads to a variation of the theoretical curves which is of the same order 



as the statistical errors on the data. The errors shown are the purely statistical binomial errors 
which are highly correlated, because all 2+1 jet events at a given y cut are included in the points 
at smaller y CVLt . The data agree with the theoretical calculations at the 15% level, however, a 
deviation from the QCD models at y cut values below 0.04 is evident in Fig. |9]. This deviation 



is correlated with the excess for z < 0.1 discussed in section O 



8 Conclusions 

The production of 2+1 jet events as defined by the JADE algorithm has been studied in deep- 
inelastic neutral current events at HERA for 160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 and 
0.04 < y < 0.95. In this kinematic range, prominent jet structures have been observed. The 
transverse energy flows with respect to the jet direction are well described by the LEPTO 6.1 
Monte Carlo program. Various measurements of the kinematic properties of the jets have been 
analysed. They are generally well described by the Monte Carlo models. The invariant mass 
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and the transverse momentum of the two jets are large enough to allow a description in terms 
of perturbative QCD. For the first time in DIS, the partonic scaling variables x p and z have 
been reconstructed from the jets and are shown to be well described by NLO calculations for 
z > 0.1. Jet rates, corrected to the parton level, have been measured as a function of y cnt 
and compared to NLO calculations. In addition to the structure function parameterisation, the 
QCD calculations have one free parameter, a s , which was taken from LEP measurements. The 
sensitivity on the choice of parameterisation is small in this kinematic regime. The dynamics of 
jet production in DIS are satisfactorily described at the 15% level by the calculations. We have 
refrained at this stage from extracting a s from our data because of the observed sensitivity at 
very small values of z and the as yet incomplete understanding of this region. 
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Figure 1: Diagrams for neutral current deep-inelastic scattering: (a) Born term, 
(b) QCD Compton scattering, (c) Boson-Gluon-Fusion leading to events with 1+1, 2+1 and 
2+1 jets, vrespectively. 
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Figure 2: Transverse energy flows in the HERA frame measured relative to the jet directions in 
1 + 1 (a-b) and 2+1 (c-f) jet events for two (x,Q 2 ) intervals (10 < Q 2 < 20 GeV 2 , 0.0012 < 
x < 0.0024) and (160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1). The data are uncorrected and are 
indicated by the full and open dots, respectively. 
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Figure 3: Transverse energy flows in the HERA frame measured relative to the jet directions in 
1 + 1 (a-b) and 2+1 (c-f) jet events for 160 < Q 2 < 1280 GeV 2 and 0.01 < x < 0.1. The data 
are shown by the full dots. The contribution from the jet which defines Arj = or A0 = for 
the data is indicated by the shaded histogram. Uncorrected data are compared to the MEPS 
and the ME simulations. 
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Figure 4: Pseudorapidity r^t of the two jets in the kinematic ranges 160 < Q 2 < 1280 GeV 2 , 
0.01 < x < 0.1 and 0.04 < y < 0.95: (a) higher 77 jet; (b) lower 77 jet. Uncorrected data are 
compared to the MEPS and the ME simulations. The boundaries of the different calorimeter 
parts are indicated. 
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Figure 5: (a) The difference in azimuthal angle, A0j et , and (b) pseudorapidity, At]^, between 
the two jets of 2+1 jet events in the kinematic ranges 160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 
and 0.04 < y < 0.95. Uncorrected data are compared to the MEPS and the ME simulations. 
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Figure 6: (a) The two-jet invariant mass V§ and (b) the jet transverse momentum Pt in the 
7*-parton centre of mass system in the kinematic ranges 160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 
and 0.04 < y < 0.95. The uncorrected data are compared to the MEPS and the ME simulations. 
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Figure 7: The distribution of z versus x v for uncorrected data in the high (x, Q 2 ) interval. The 
area defined by the curve indicates the kinematic limit for y cut = 0.02 in the range 0.01 < x < 
0.1. The limit on z is given by z min ~ y cut = 1 - z max . 
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Figure 8: (a) Uncorrected x p distribution of two-jet events and (b) uncorrected z distribution 
of the jet with the smaller z. The uncorrected data are compared to the MEPS and the 
ME simulations. The shaded histogram in (a) shows the uncorrected x p distribution in the 
kinematic ranges 10 < Q 2 < 20 GeV 2 , 0.0012 < x < 0.0024 and 0.04 < y < 0.95. 
(c) corrected x p distribution and (d) corrected z distribution compared to the NLO calculation 
(PRO JET). The data have been corrected with the MEPS simulations. Statistical errors only 
are shown. 
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Figure 9: The corrected jet production rate R n +i in % as a function of y cnt is compared to 



the NLO calculation of the programs PROJET [IT) and DISJET ]L6] in the kinematic ranges 
160 < Q 2 < 1280 GeV 2 , 0.01 < x < 0.1 and 0.04 < y < 0.95. The data have been corrected to 
the partonic level with the MEPS model. Statistical errors only are shown. 
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